Helicity order: hidden order parameter in URu2Si2 
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We propose that the "hidden order parameter" in URu2Si2 is a hehcity order which must arise, if 
the Pomeranchuk criteria for the spin-antisymmetric Landau parameters with respect to the stabihty 
of a Fermi hquid state are violated. In a simple model, we calculate the specific heat, linear and 
nonlinear magnetic susceptibilities and the change of transition temperature in a magnetic field 
with such an order parameter, and obtain quantitative agreement with experiments in terms of two 
parameters extracted from the data. The peculiar temperature dependence of the NMR linewidth 
and the nature of the loss of excitations in the ordered phase seen by neutron scattering are also 
explained and experiments are suggested to directly confirm the proposed order parameter. 

PACS numbers: 71.27.+a, 71.10.Hf, 75.10.-b, 75.40.Cx 



The "hidden order" phase in the heavy fermion com- 
pound URu2Si2 below the second order transition at 
17. 5K l] has remained a puzzle for about 20 years. The 
magnitude of the specific heat at the transition is equiv- 
alent to that of ordering of a moment of about 0.5//_b per 
unit cell. No change in spin-rotational symmetry or lat- 
tice translational symmetry consistent with this specific 
heat has been discovered. Detailed neutron diffraction 
experiments^ reveal a moment of only about 0.03/i 
per unit cell, which as NMR and /iSR experiments 
reveal, is due to the presence of a second phase. Some 
very interesting proposals for new types of order have 
been made 0, |^,1Bl] i which have not been supported by 
experiments designed to look for them. 

Some of the other properties measured at the transi- 
tion to the "hidden order" phase and in it are, the linear 
magnetic susceptibility which only changes slope at the 
transition, the nonlinear magnetic susceptibility which 
shows a singularity at the transition similar to that of 
the specific heat the change of the transition tem- 
perature with an applied magnetic field H, , the loss 
of low energy excitations observed by neutron scatter- 
ing for a range of wave vectors |2, llfl 113 and in trans- 
port measurements'!^, and the NMR relaxation ratefisj 
which exhibits the extraordinary result that there is an 
extra inhomogeneous relaxation rate in the ordered phase 
which increases below the transition temperature propor- 
tional to an order parameter. 

We suggest here that the transition is to a state 
proposedjij as a cure to the spin-antisymmetric Landau- 
Pomeranchuk instability(LPI) of the Fermi-liquid. For 
reasons which will be clear, we call such states helicity or- 
dered states^^. We calculate the thermodynamic prop- 
erties near the transition, and account quantitatively for 
the observed thermodynamic features and qualitatively 
for the NMR and the excitation spectra with parameters 
extracted from the experiments. We also suggest experi- 
ments which can provide direct evidence for the proposed 
phase. 

In Landau's Fermi liquid theory the change in the 
free-energy due to a a small change of the equilibrium 



distribution function (5n(kcr) is 
5F = ^e^(5n(ko-) 

kcr 

- /(kcr, k'cr')(5n(kCT)(5n(k'CT') 
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The interaction functional /(k(T,k' cr') has 
symmetric(s) and spin- antisymmetric (a) parts: 



spin- 



/(ka, kV) = /^(ka, kV) + /«(ka, k'a'). (2) 

The coefficients of an expansion of f^'"^ in terms 
of the irreducible representations of the Fermi-surface 
are the Landau parameters F^'"" , in terms of which 
Pomeranchuk [T^ obtained a set of conditions for the 
stability of the Fermi-liquid: I + [21 + > 0. 

Any violation of these conditions leads to a Landau- 
Pomeranchuk instability, which must be cured by a bro- 
ken symmetry in corresponding irreducible representa- 
tion / and spin symmetry s or a. For example, the 
ferromagnetic instability occurs for Fq —>■ —1. Spin- 
symmetric instabilities in finite /-channels have attracted 
much recent interest^J^ 0, 0|. A spin-ordered state, 
which is anisotropic in momentum space(without change 
in translational symmetry) is the obvious cure to the fi- 
nite I antisymmetric LPIs [T3. l20l |. We develop this idea 
here; we find that besides the LPI criteria, additional 
condition must be satisfied so that the instability is of 
second order. 

Consider the model Hamiltonian, 



H 



ka 

\ Jk,k'(q) (cl+gO'Cfc) • (4 
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where el is the spectrum of a nonintcracting Fermi gas, a 
Pauli matrices, q is the momentum transfer; of interest 
is the instability in the forward scattering limit ((7 — > 0). 
Jk,k' is the interaction in spin- antisymmetric channels, 
which can be expanded as Jk,k'(0) = J2i •^/-Pi(cos0k,k'): 
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FIG. 1; Schematic Fermi surface for the two different spin- 
directions. The dotted line denotes the Fermi surface of the 
paramagnetic phase with vanishing order parameter while the 
solid lines illustrate the Fermi surface of the helicity ordered 
phase. 



where 9u,k' the angle between k and k', and Pi {x) Leg- 
endre polynomials. 

In the normal state of a Fermi-liquid, helicity is disor- 
dered since the spin-quantization axes at each k can be 
independently rotated. The proposed order parameter 
for the model has the general formjT3.l2(Tl|: 



((5n(k,cr)) =cr •D(k/). 



(4) 



The spin-quantization axis is thereby fixed in relation to 
the direction on the Fermi surface |2ll|. 

We need consider only one specific / channel and write 
Jk.k' in a separable form J/ P/ (cos (cos ^fc'). The sim- 
plest order parameter has 13^ (k/) ^ 0, so that the asso- 
ciated energy parameter is 



Ai = {JiaW{kf)) = C^JiPiicosek)ink^-nki)). 

I. 



(5) 



This Ising order parameter is especially useful to discuss 
the tetragonal compound URu2Si2, which has a large 
anisotropy in the magnetic susceptibility favoring the c- 
axis. With Eq. @ we have a noninteracting model with 
the effective spectrum 



(6) 



where H is the external magnetic field, and no the ef- 
fective single-electron magnetic moment. Imposing the 
requirement of a constant chemical potential, the Fermi 
surface for the up and down spins are split as schemati- 
cally illustration in Fig^ 

We calculate the free energy following standard meth- 
ods. The following results are for / = 1, but can be easily 
generalized to higher-Z channels. The free energy can be 
in general separated into two parts, 



(7) 



where describes a paramagnetic phase(Ai = 0). 
The specific heat and the magnetic susceptibility [M — 
XiH+{x3/-^^-)H'^ + . ■ ., Xi and X3 are the linear, nonlinear 
spin susceptibility, respectively] can be easily obtained. 
Including terms of 0{H^), 0{T'^) and the variation of 



the density of states p{e) near the chemical potential to 
0[p"{eF)\, we obtain 
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When H = and T = 0, the criterion to have a con- 
tinuous phase transition is A < and B > 0, i.e., 

„'2 
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This, together with Eo. lfTmi . requires p" / p > p'^ / p^ ■ 
presence of a small magnetic field, varies as 



H 



T, 



Ho 



Te [1 - [H/Hof] , 
(2/3)i/27rfcBTe/Mo 



1/2 ■ 
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For noninteractingelectrons, these are standard re- 
sults(see, e.g., Ref.(2a|'): we list them here to use them 
to extract parameters from the normal state experimen- 
tal results. For interacting electrons (in the limit of zero 
field) they are multiplied by Landau parameters, m*/m 
for the specific heat, [ra* /m) / [l+F^) for the susceptibil- 
ity and (m*/m)/(l+4FQ ) for the nonlinear susceptibility. 

Tm{T, H,Ai), is the additional contribution to the free 
energy for Ai 7^ 0. Expressing it in series of the order 
parameter Ai gives. 



(9) 



(10) 



The first gives Ji < and p\Ji\ > (2/ + l)/2 (here I = 1), 
which is precisely the LPI criterion. The second is an 
additional criterion, on the form of the density of states 
at the Fermi surface to have a second order transition. 
li B < 0, one must expand the free energy to terms of 
order Af. In that first-order phase transition is 

favored. 

At iJ = 0, the critical temperature is given by 



--^ • (11) 



In 



(12) 
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where h, ^ [p"/p - p'^ {3p^)]/ {p" / p - p'^p^). 

Below Tc we have the non-trivial solution = —A/B, 
i.e., a helicity-ordered state. The changes in some ther- 
modynamical quantities at and below Tc from their val- 
ues for T > Tr are calculated to be 



Sch/c°h 



9 xr 

2p pI 



q9p 

X3 



(13a) 
(13b) 
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where g, ^ {p"/p - p'^ p^ /[p" /(bp) - p'V(3p^)]. The 
specific heat shows of-course a characteristic mean-field 
discontinuity at the transition point; more interesting is 
the fact that X3 also shows a discontinuity while the linear 
magnetic susceptibility shows merely a change of slope. 
A singularity in the nonlinear magnetic susceptibility is 
to be expected for any order parameter O when a term 
|0pi?2 is allowed in the free energy. Usually the coeffi- 
cient of this term is so small that the singularity in X3 
is not noticed. What is special about URu2Si2 is that 
the dimensionless mean-field jump in X3 is similar to the 
dimensionless mean-field jump in the specific heat. This 
and and several other properties are quantitatively ex- 
plained below. 

We have ignored the Landau parameters in the dimen- 
sionless quantities in Eas.(|51 [T^ . The reason is that the 
change in the Landau parameters near the transition may 
be shown following Leggett [2^ , for the case of transition 
in superfluid ^He, to be proportional to (pAi)^. To this 
order, assumed ^ 1, they vanish in dimensionless quan- 
tities. 

In the following, we try to fit the experimental data: 
the specific heat from Ref. , and the linear and nonlin- 
ear magnetic susceptibilities from Ref.Q- The data are 
shown in Fig|21for the reader's convenience. From Eqs. 
^ and (|13|l . in addition to the prefactors, 70 = ^^pkg, 
Xo — 2pgp, xo — 3!/iQP, all other quantities can be deter- 
mined by two additional dimensionless variables, Ci = 
{p"/p-p'yp^)iny6)ikBT,)^ and C2 - (p'Vp')/(p7p), 
where Tc is taken as 17.5K. 

Consider the linear magnetic susceptibility. It is con- 
tinuous with a change in the slope at the transition point, 
which is consistent with the result in Ea. ljl3b(l . Also no- 
tice that in the normal state it shows a significant linear 
temperature dependence besides the constant Pauli term. 
This is indeed required by the theory in order that there 
be a second order transition to the helicity ordered phase 
[see Eq. (Unj]. Near 

Xi{T>Tc) = xo(l + Ci) + 2xoCi(T-r,)/T„ 
SxiiT<Tc) = xoyCi^l^^(r-T,)/r,. (14) 

By fitting them with the experimental data, we can ex- 
tract 



e 

g 200 



(a) 




01 


(c) 








emu/mc 










io . 















(b) 



(d) 



FIG. 2: The experimental data of direct relevance to the cal- 
culations in this paper. The specific heat data(a) are ex- 
tracted from Fig. 1 in Ref.^. In the specific heat data a con- 
tribution to C/T proportional to T^, presumably mostly due 
to phonons, has been subtracted off. Data of the linear(b) and 
nonlinear(c) magnetic susceptibilities (along the c-axis) are 
extracted from Fig. 2 in Ref. [3 • (d) shows the inhomogeneous 
linewidth of Si-NMR for magnetic fields in the c-direction and 
in the plane, which is extracted from Fig. 4 in Ref. 13]; the 
solid line is the fitting function A = 12[1 - {T/TcfY'/'^{G). 



(with Xo ~ 63.0 emu/mole T). We then can calculate 
other thermodynamic quantities. At the disconti- 
nuities of the specific heat coe0icient(7 — C/T) and 
the nonlinear magnetic susceptibility are calculated from 
Eq.lO to be 



5l{Tc)ho = 1.4, 5x3{Tc)/x3{T^c 



2.1, 



(16) 



while the corresponding experimental quantities from 
Fig. 121 are approximately 1.5 and 2.4, respectively. In 
finite magnetic field, Ea. (|12|l predicts that Tc decreases 
as (H/Hq)'^, where Hq — 38.2T if we take po to be one 
Bohr magneton. This again, is in agreement with the ex- 
periments, where Hnis estimated as 48.5(1)T in Ref.|l], 
while 35. 3T in Ref.j^]- We can also determine the order 
parameter {H = 0) 



Ai(r) - 77 [1 - [T/Tcf] 



1/2 



K. 



(17) 



Ci « 0.35, C2 « 0.52, 



(15) 



To summarize, we get the qualitatively correct behav- 
ior of the linear susceptibility and extracting two parame- 
ters from it, can explain quantitatively the relative jump 
in the specific heat and that of the nonlinear suscepti- 
bility as well as the characteristic field for the suppres- 
sion of the transition. (The fact that a simple model for 
the Fermi surface of the paramagnetic phase gives these 
quantities quite well is doubtless due to the fact that we 
are comparing dimensionless quantities.) We do not do 
well on the specific heat just below Tc] our slope is about 
a factor of 2 smaller than the straight line fit one might 
make in Fig|21 Note that the proposed order parameter 
has an Ising symmetry. So, the actual exponents in the 
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critical regime of the transition are expected to be that of 
an Ising model in 3-dimensions; for example the specific 
heat is expected to show a lambda shape. While further 
experiments are required to test this, the measured spe- 
cific heat is not inconsistent with such a form (see Fig. 
EJ. In such a case, a mean-field fit to the data always 
gives a slope smaller than the experiment. 

The values of Ci, C2 required above imply that p" / p ~ 
2p'^/p^, i.e., the chemical potential in the normal phase 
lies near a local minima of the density of states. This is 
consistent with the density of states calculated by band- 
structure calculations p^. 

Let us next consider the NMR measurements. Si- 
NMR[31 reveals no change in Knight shift but an in- 
crease in the inhomogeneous linewidth below Tc which 
within experimental uncertainty can be fitted to be cx 
[1 - (T/rc)2]i/2[see FigOd)], i.e. proportional to an or- 
der parameter. This is quite unusual. We note first that 
in a perfectly pure sample, we expect no change in Knight 
shift or linewidth. In the presence of impurities, which 
locally break the reflection symmetry about the basal 
plane, a local ferromagnetic region forms, as is evident 
from FigQ] The magnitude of the local fleld then is pro- 
portional to the order parameter but its magnitude as 
well as direction are random. This gives no Knight shift 
but a linewidth consistent with observations. The mag- 
nitude depends on details of the defect, but a O(0.5/iB) 
defect, expected from the magnitude of the order param- 



eter, need be present only in concentrations of a few parts 
in a thousand to produce the observed linewidth of order 
10 Gauss. The observed linewidth is almost independent 
of the direction of the applied magnetic field. This can 
be shown to occur for generic distribution of impurities 
about the Si-sites psj. 

We intend to calculate the excitation spectra in the fu- 
ture. One can however see that, given the Fermi-surfaces 
shown in FigQ] a decrease in inelastic scattering below 
a certain characteristic energy of the order of Ai is ex- 
pected for spin-fiip particle-hole excitations. This is what 
is observed in neutron scattering Ullillll. A quanti- 
tative calculation requires a more realistic model of the 
normal state Fermi-surface than used here. We note that 
the characteristic magnitude of this energy scale observed 
in inelastic neutron scattering experiments is0 ~ TO-fC, 
which compares well with our estimate [see the coefficient 
in Eq.ini)]. 

Finally, we turn to how the proposed order parameter 
may be directly observed. On applying an electric field 
parallel to the c-axis, a spin-current would be generated 
but no such effect should occur on applying electric field 
parallel to the basal plane pTj. Another direct possibility 
is thr oug h spin-polarized positron annihilation suggested 
to uspq. 

We wish to acknowledge useful discussions with V. Aji, 
W. Buyers, D. MacLaughfin, A. Mills, A. Mighori, J. 
Mydosh, and A. Ramirez. 
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